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At numerical resolutions around 512 3 and above, three-dimensional energy spectra from tur- 
bulence simulations begin to show noticeably shallower spectra than k~ 5 ^ 3 near the dissipation 
wavenumber ('bottleneck effect'). This effect is shown to be significantly weaker in one-dimensional 
spectra such as those obtained in wind tunnel turbulence. The difference can be understood in 
terms of the transformation between one-dimensional and three-dimensional energy spectra under 
the assumption that the turbulent velocity field is isotropic. Transversal and longitudinal energy 
spectra are similar and can both accurately be computed from the full three-dimensional spectra. 
Second-order structure functions are less susceptible to the bottleneck effect and may be better 
suited for inferring the scaling exponent from numerical simulation data. 
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I. INTRODUCTION 

Based on dimensional analysis, Kolmogorov con- 
cluded that the energy spectrum for isotropic hydrody- 
namic turbulence has the form 

E(k) = C K e 2/3 £T 5/3 (1) 

for wave numbers k in the inertial range between the 
energy-carrying and the dissipation wave number fed = 
e 1 / 4 ;/ -3 / 4 (where v denotes the kinematic viscosity, e the 
spectral energy flux, and Ck is nowadays called 'Kol- 
mogorov constant'). This scaling has been confirmed ex- 
perimentally over several orders of magnitude [2>kj- Nev- 
ertheless, numerical simulations consistently show excess 
power just before the dissipation wavenumber k<u which 
manifests itself particularly at high resolution Q. This 
phenomenon has been named 'bottleneck effect' 0, Q 
and is usually explained by the lack of smaller-scale vor- 
tices at wave numbers k > k^, which makes the energy 
cascade less efficient around fed. According to a related 
interpretation, the bottleneck effect is the consequence of 
viscosity stabilizing small vortex tubes against the kink 
instability 0. 

The effect is particularly strong when an unphysical 
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hyperviscosity is used; see Refs [alSj for results from two- 
dimensional hydromagnetic turbulence. In experimental 
data, on the other hand, the bottleneck effect is less pro- 
nounced 0, |n| , and it has previously been noticed that 
it is much weaker when one-dimensional energy spectra 
are used . In the present paper we discuss a simple rela- 
tion between one-dimensional and three-dimensional en- 
ergy spectra, which agrees well with the simulation data 
and explains this difference. 

The data we use for discussing the bottleneck effect are 
from a weakly compressible isothermal three-dimensional 
forced turbulence simulation at a numerical resolution of 
1024 3 grid points. The forcing has vanishing net helicity 
and the forcing wave number kf is between 1 and 2. The 
box size is L x = L y = L z = 2n, which discretizes the 
wave numbers in units of k\ = 1 . The viscosity v is cho- 
sen such that the Reynolds number based on the inverse 
mean forcing wave number, u rms / \vkf), is around 1700. 
The Taylor microscale is \/5w rms /o; ri ns ~ 0.14, where 
u) Tms is the root mean square vorticity, so the correspond- 
ing Taylor microscale Reynolds number is 350. The av- 
erage dissipation rate e is such that fcd/fcf is around 130. 
The root mean square Mach number is between 0.17 and 
0.20; for this type of weakly compressible simulations, we 
find that the energies of solenoidal and potential compo- 
nents of the flow have a ratio E pot /E so \ 10~ 4 -10~ 2 for 
most scales; only towards the Nyquist frequency the ratio 
increases to about 0.1. Even for Mach numbers between 
0.5 and 10, this ratio is only about 0.1 to 0.2 [T2I Il3|. 
It is thus reasonable to assume that compressibility is 
irrelevant for the results presented here. This is also sup- 
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FIG. 1: Comparison of the averaged one-dimensional lon- 
gitudinal and transversal spectra, Ej J (k) and Er(k), re- 
spectively, with the three-dimensional spectrum E(k) for a 
forced turbulence simulation at 1024 3 grid points. Top: 
the spectra. Bottom: 'compensated spectra' e~ 2 ' 3 k 5 ' 3 E{k), 
e~ 2/3 k 5/3 E L/T (k); the horizontal line represents a k~ 5/3 Kol- 
mogorov spectrum. The local maximum of E(k) around 
k ps 30 represents the bottleneck effect. The dissipation wave 
number is k d = e 1 / A v~ 3 / A w 200. 



one-dimensional spectra have been averaged horizontally, 
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where |u±| 2 = Iw^j 2 + \u y \ 2 . Here 



u(x,y,k) 



LA 1/2 1 N ' 
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J2e lkZn u{x,y,z n ) (4) 



n=l 



is the one-dimensional Fourier transform of the veloc- 
ity vector u(x), and k scans the interval [0, &N y ] = 
[0,tt/8z] at the resolution 5k — 2n/L z . The three- 
dimensional spectrum E(k) has been obtained by in- 
tegrating the three-dimensional spectral energy density 
over shells k% — 5k/ 2 < k < ki + 5k/ 2 in analogy to 
Eq. (|A~fi)) . with 5k = 2tt/L x (and L x = L y = L z ). The 
spectra are normalized such that 

f°° 1 
jf E(k)&k = -(u 2 ) 

= 3 E h (k)dk = 3 / E T (k)dk. (5) 
Jo Jo 



II. RELATIONS BETWEEN ONE- AND 
THREE-DIMENSIONAL SPECTRA 

While most experimental measurements yield lon- 
gitudinal one-dimensional spectra E^(k), the discus- 
sion of the relation between one-dimensional and three- 
dimensional spectra is significantly simpler for the total 
one-dimensional spectrum E\Y>{k). We thus split this 
section in two parts: we first outline the relations for 
the total one-dimensional spectrum Eyd (k) , while in the 
second part we obtain analogous results for longitudinal 
spectra, which are of direct relevance for experiments. 



ported by the fact that incompressible pseudo-spectral 
simulations at a resolution of 1024 3 show a bottleneck 
effect |j . 

The simulations discussed here were carried out us- 
ing a high-order finite-difference code [l4| and thus 
complement the results that have so far been ob- 
tained using spectral codes. Figure Q shows the three- 
dimensional spectrum E{k) together with the longitudi- 
nal and transversal one-dimensional spectra Ei,{k) and 
E^{k) in the upper panel. The 'compensated spectra' in 
the lower panel allow easy identification of the bottleneck 
effect and show that it is practically present only in the 
three-dimensional spectrum E(k). 

In order to reduce the otherwise huge fluctuations, the 



A. The total one-dimensional spectrum 

The total one-dimensional spectrum En)(k) is the sum 
of the longitudinal and twice (for the two directions) the 
transversal one-dimensional spectra, 

E 1D (k) = E L (k) + 2E T (k). (6) 

It is thus in some sense 'more isotropic' than its con- 
stituents which results in simpler relations to the fully 
isotropic three-dimensional spectrum. In Appendix IT1 we 
show that Em(k) is related to E(k) by 

E 1D (k)=f k °°^dk\ (7) 
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which can also be inverted to give 



and 



E(k) 



dE 



ID 



dk 



(8) 



provided the turbulence is isotropic. As E(k) must be 
positive, Eq. J7J) shows directly that the one-dimensional 
spectrum Eny(k) must be monotonously decreasing. On 
the other hand, no such restriction holds for E(k), which 
is in fact increasing near k — 0, since E(0) = 0. 

Equation JHJ is a local relationship, and thus the func- 
tional form of E(k) is fully determined by the local be- 
havior of Em(k) at a given wave number. To relate it 
to the bottleneck effect, we introduce the compensated 
spectra 

£ 1D (fc) = £; 5 / 3 £ 1D (fc), E{k) = k 5 / 3 E(k), (9) 
and rewrite Eq. (jSJ) in the form 



E(k) 



dln£a D (fc) 
dlnfc 



E m (k). 



(10) 



Several conclusions can be drawn from Eq. (|10fl . First, if 
there is a finite interval where the Kolmogorov scaling 
holds for the one-dimensional spectrum, then the same 
scaling will hold for the three-dimensional spectrum, and 
vice versa. Further, Eq. 1)10(1 shows explicitly that even 
if EiY)(k) is monotonously decreasing, E(k) may show a 
maximum near fed, provided that Em(k) bends towards 
the dissipative range sufficiently suddenly. For example, 
if the one-dimensional spectrum has the form 

E m (k) = fc- 5 / 3 exp[-(fc/fc d )"], (11) 

then for small wave numbers k the compensated three- 
dimensional spectrum behaves like 



E(k) 



(for k < fc d ), (12) 



5 / 5\ / k 

which shows that there will be a bottleneck effect in the 
three-dimensional spectrum when the one-dimensional 
spectrum falls off with an exponent n > 5/3. Finally, 
if Em(k) shows a bottleneck effect, i.e. En)(k) shows a 
local maximum at some wave number fc m , then Eq. I|10|) 
shows that E(k) will also have an enhanced value there, 
E(k m ) = (5/3)E 1I) (k m ), and thus the three-dimensional 
spectrum E(k) must show a bottleneck effect, too. 



B. Longitudinal and transversal one-dimensional 
spectra 

For isotropic turbulence the one-dimensional longitu- 
dinal and transversal spectra are related to the three- 
dimensional spectrum by (see Appendix ^J) 



i r°° / k 2 



E(k') 

k' 

E(k') 



dk' 



dk' 



(13) 
(14) 



E{k) = k 2 E'l(k) - kEl(k) 



(the primes denoting derivatives). 
Eq. i|13|) . we obtain 



(15) 
Differentiating 



EL(k) = -k J°°^p-dk'<0, 



(16) 



which shows that E^k) (just like Exoik)) must be 
monotonously decreasing. 

Equation H15|l is again a local relationship, and thus 
the functional form of E(k) is fully determined by the 
local behavior of E\^{k) at a given wave number. As in 
Eq. jnj, we introduce the compensated spectra 

E(k) = k 5 / 3 E(k), E h (k) = k 5 / 3 E L (k), (17) 



and rewrite Eq. (|15J) in the form 
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55 
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E, 



13 
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k 2 El 



(18) 



The conclusions that can be drawn from Eq. (|18|) are sim- 
ilar to those for En)(k). If there is a finite interval where 
the Kolmogorov scaling E^(k) = C^e 2//3 fc _5 / 3 holds for 
the longitudinal spectrum (which implies E^ = E^ = 0) 
then the same scaling will hold for the three-dimensional 
spectrum, and vice versa. The Kolmogorov constant 
Ck of the three-dimensional spectrum is then a factor 
55/9 « 6.1 larger than the corresponding constant 
for E h {k)- cf. Eq. (16.37') of Ref. jf^. 

Equation (|18l) shows that even if E-^(k) is 
monotonously non-increasing, E(k) may show a 
maximum near fed, provided that the term involving 
E^ dominates over the curvature term k 2 E^. Finally, 
if E^(k) shows a bottleneck effect, then beyond the 
local maximum, E[^(k) < 0. If the maximum is wide 
enough, Eq. 1)18(1 implies that E will show a bottleneck 
effect, too. The situation is less clear if k 2 a E^(k m ) is 
large (corresponding to a narrow maximum of E^(k)), 
but numerical experiments with model spectra suggest 
that the three-dimensional spectrum E^(k) will still 
be non-monotonous, although it may vary in a more 
complicated manner. 



III. SAMPLE SPECTRA 

A. Model spectra 

For illustration purposes, we consider a longitudinal 
spectrum of the form 



-5/3 

E h (k) =[-£-) exp[-(k/k d y 



(19) 
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FIG. 2: Compensated three-dimensional spectrum E(k) = 
(fc/fcd) 5//3 E(k) corresponding to the longitudinal spectrum 
(1191 1 for exponents n — 1,2,4. The bottleneck effect appears 
for n > 2. 



The compensated spectrum E-^(k) = (k/kd) 5 ^ 3 Ei,(k) = 
exp[— (fc/fcd)"] is monotonously decreasing and thus 
E-^(k) shows no bottleneck effect at all. At k = fed the 
compensated three-dimensional spectrum has the value 



E(k d 



55 



1G 



(20) 



which shows that E{k d ) > 55/9 (a sufficient condition for 
the bottleneck effect in the three-dimensional spectrum) 
if n > (55/48)(e - 1) w 1.97; this is illustrated in Fig. 
A more thorough analysis reveals that E(k) will have a 
maximum if, and only if n > 5/3, but for n < 2 the 
maximum is hardly discernible. 

As we have seen, in order to get a bottleneck effect, we 
need a one-dimensional spectrum with a more compli- 
cated form than just fc~ 5 / 3 exp(— fc/fcd). One functional 
form where this is given has been proposed by She and 
Jackson [Tfjj |. based on experimental data: 



E L (k) 
Ej J (k p ) 



-5/3 



+ 0.8 



3 -0.63fc/fc p 



Figure Efa) shows this spectrum [for Ez,(k p ) = k v 



(21) 
= 1], 



together with the corresponding three-dimensional spec- 
trum (|15|1 ■ The compensated spectra clearly show a bot- 
tleneck effect in the three-dimensional spectrum, which 
is (practically) absent in the one-dimensional spectrum. 

Quian has proposed a spectrum that shows quite 
a marked bottleneck effect. Based on a closure model, he 
suggests the functional form 



E(k) 



£ 2/3 fc -5/3 



1.19 + 6.31 



2/3* 



x exp[-5.4(fc//fc d ) 4 / 3 ]. 



(22) 
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FIG. 3: Comparison of compensated spectra Ei,(k) and E(k) 
for the models of She & Jackson, and of Quian. (a) longitu- 
dinal spectrum l!2H , together with derived three-dimensional 
spectra, longitudinal and total ID spectra; the spectra are 
normalized according to k p = E^(kp) = 1. Note the appear- 
ance of a mild bottleneck effect, i.e. a maximum in E(k). (b) 
Three-dimensional spectrum I22H . together with derived one- 
dimensional spectra; the spectra are normalized according to 
fed = e = 1. The bottleneck effect is quite pronounced and 
appears in both, one- and three-dimensional spectra. Note 
that in both plots the amplitudes of the one-dimensional spec- 
tra have been scaled to get matching plateaus in the inertial 
range. The dotted horizontal lines are drawn for orientation. 



Figure[3Jb) shows the longitudinal and three-dimensional 
spectra. For this model, the bottleneck effect is very 
prominent in E(k), and also evident (although weaker) 
in the longitudinal one-dimensional spectrum E-^{k). 
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B. Spectra from direct numerical simulations 

1. Three-dimensional spectrum from Ej J (k) 

To avoid (double) numerical differentiation of our spec- 
tra, we use a parameterization for E\ J (k). The one- 
dimensional spectrum shown in Fig. ^ is weu approxi- 
mated by the formula 



E^(k) = (fc- 5 / 3 + a fc ai K Q2 * ° 3 



+ 0,4k + a^k 2 
1 + a 6 k + a-jk 2 



(23) 



with Oi = (2.1 xl(T 8 , 4.3, 0.42, 0.85, 1.2,0.00048, 
—0.0068,0.34), where k = k/k p . The peak dissipation 
wavenumber k p « 18 is the location of the maximum of 
the dissipation spectrum k 2 E^^ and is about one order of 
magnitude smaller than fed |H|. Figure 0] shows the pa- 
rameterized one-dimensional spectrum E^(k) according 
to Eq. 112311 . together with the derived three-dimensional 
spectrum E(k) from Eq. i|15|) . Also shown are data 
points from the numerical simulation (diamonds and 
crosses). Comparing the calculated three-dimensional 
profile (solid line) with the data points (crosses), we 
find that Eq. (|15|l agrees quite well with the numerical 
data for not too small k. The discrepancy for very small 
wave numbers can be explained by the fact that the 
periodicity of the numerical box precludes isotropy at 
the largest scales. 

It is quite evident from Fig.0]that E{k) shows a much 
more pronounced local maximum near the dissipation 
wavenumber k^ than does E^(k). The width and struc- 
ture of the maximum is well-reproduced by Eq. (|15fl ap- 
plied to our parameterization (|23|l of E-^(k). 



2. One-dimensional spectra from E(k) 

In contrast to experiments, data from numerical sim- 
ulations easily provide the three-dimensional spectrum 
E(k). Sometimes it may be interesting to determine 
the one-dimensional spectra E^,(k) and £r(fc) from the 
three-dimensional spectrum. The corresponding rela- 
tions lO, en were given above. Here we apply them 
to our numerical data to see how well the inferred one- 
dimensional spectra agree with those directly obtained 
from the simulation data. Figure |3] shows the longi- 
tudinal and transversal spectra obtained by applying 
Eqs. i|13[l and i|14|) to the three-dimensional spectrum us- 
ing the trapezoidal rule. The agreement is quite good for 
both one-dimensional spectra apart from the very lowest 
wave number. 
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FIG. 4: Compensated longitudinal energy spectrum 
k 5 ^ 3 E^(k) according to the parameterization 123t and corre- 
sponding three-dimensional spectrum k 5 ' 3 E(k) from Eq. 1150 . 
For comparison, crosses and diamonds show the spectra ob- 
tained directly from the simulation. The spectra have been 
normalized by introducing k = k/kp, El = El I Eh(kp) and 
E = E I Ei J (k p ) , where k p is the peak dissipation wavenumber. 
The three-dimensional energy spectrum agrees quite well with 
the numerical data for k > 45k, where Sk = 2n/L x ~ 0.06 k p 
is the wave number resolution. 
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FIG. 5: Compensated one-dimensional energy spectra. 
Shown are the longitudinal spectrum E^(k) and the transver- 
sal spectrum ET(k). Data points obtained directly from the 
numerical simulation are indicated as boxes (El) and crosses 
(Ex)- The dashed and dotted lines show the one-dimensional 
spectra obtained from Eqs. II -'ill and (1141 . respectively. The 
spectra agree quite well with the numerical data for k > 2 5k. 



IV. STRUCTURE FUNCTIONS structure function, 

Another classical tool in turbulence research is the in- 
vestigation of structure functions @. The second-order ^(r) — (i u ( x ) — u(x+r)] 2 ") , (24) 
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is related to the three- and one-dimensional spectra via 
the Fourier-type integral transforms [cf. fl5| 



(25) 



= 4 y (3 — 3 cos kr + kr sin kr) El (k) dk . (26) 
o 

Similarly, the longitudinal and transversal second-order 
structure functions 



Si h) (r) = {[u z {x)-u z {x+rz)f), 
si T) {r) = ([u z (x)-u z (x+rx)} 2 ), 
can be expressed as 

oo 

S< L/T) (r)=4 J (1 - cos kr)E L/T (k) dk. 
o 



(27) 
(28) 



(29) 



Localized variations of E(k) or E^k) in wavenumber 
space, such as the bottleneck effect, will influence S^M 
in a nonlocal fashion. Correspondingly, lit- 
tle insight into the bottleneck effect can be expected from 
structure functions. On the other hand, this means that 
structure functions are less sensitive to the bottleneck 
effect and might thus be a more robust tool for assess- 
ing scaling exponents and possibly even the Kolmogorov 
constant at moderate Reynolds number. 

However, while structure functions are much smoother 
than spectra, their scaling range is considerably smaller, 
which adds its own difficulties to that method. In 
three numerical simulations at 256 3 , 512 3 and 1024 3 grid 
points, we find the value of the structure function ex- 
ponent derived from ^(r) to be 0.74, 0.67, and 0.68, 
respectively, which is quite close to Kolmogorov's value 
of 2/3. However, if we use the the two transversal and 
the longitudinal structure function (corresponding to the 
^-displacement of the u x , u y and u z component of the 
velocity vector) , the values are far less accurate and span 
the intervals [0.67,0.79], [0.62,0.72], [0.60,0.72], respec- 
tively, for the three resolutions, which indicates that the 
convergence is maybe not yet reached. 

The situation for the Kolmogorov constant is even 
worse, because from S^M we get the values 3.96, 1.93, 
1.8, respectively, for the three resolutions. So in practice, 
the slow convergence may well render this method more 
unreliable than the use of one-dimensional spectra. 

Similar to the spectra, the second-order structure func- 
tions can be transformed into each other (assuming 
isotropy) according to the relations 



S 2 (r) 



r dSl^r± +35 (L )(r)) (30) 



dr 



Si h) (r) = 1 ^ S 2 (r')r' 2 dr', 
r Jo 



(31) 




(r)/4.82 
11/3 (£r)- 2/3 4 T> (r)/4.82 
11/4 (zry m $Xr)l 4.82 



0.01 



0.10 



1.00 



k,r 



FIG. 6: Compensated second-order structure functions from 
our direct numerical simulation data. The amplitudes have 
been scaled to get matching plateaus representing the Kol- 
mogorov constant Ck- 



sP(r) 



S 2 (r)-Si h \r) 



(32) 



Like in the case of relations l|13ll - (|I5[) . this has implica- 
tions for the monotonicity of the compensated structure 
functions S(r) = SY~ 2 / 3 (r) near the boundaries of the 
scaling interval; see also Ref. Normally, however, 

none of the structure functions S^r), 



show a secondary bump near the edges of the scaling in- 
terval, as can be seen in Fig. which implies that we 
cannot directly apply the results of Sees. Inland Sec ITTTl 
to structure functions. 



V. CONCLUSIONS 

In this paper we have highlighted the discrepancy be- 
tween one-dimensional and three-dimensional spectra. 
Well-known relations between the spectra show that the 
three-dimensional spectrum may show the bottleneck ef- 
fect even if the one-dimensional spectra do not show it 
at all, while the converse cannot happen. The spectra al- 
ways agree in the inertial range, E(k) cx Ej J (k) cx fc~ 5 / 3 , 
but in current numerical simulations the length of the in- 
ertial range is limited to about one decade, so the discrep- 
ancy is quite noticeable. Indeed, the topic of a bottleneck 
effect in hydrodynamic turbulence has only emerged in 
the last ten years since numerical simulations have shown 
this to be a strong effect. On the other hand, the relation 
linking one-dimensional to three-dimensional spectra has 
been known for fifty years, but it has to our knowledge 
never been explicitly discussed in connection with the 
bottleneck effect. In the present paper we have shown 
that much of the bottleneck effect seen in numerical tur- 
bulence simulations is simply the result of the mathe- 
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matical discrepancy between one- and three-dimensional 
spectra. 

The bottleneck effect has no evident manifestation in 
the second-order structure functions, where localized fea- 
tures in fc-space appear in a delocalized manner. This 
implies that in order to obtain the asymptotic energy 
spectrum exponent it may be easier to use second-order 
structure function exponents, although in practice the 
reduced scaling range may render this method difficult. 



APPENDIX: RELATIONS BETWEEN ONE- AND 
THREE-DIMENSIONAL SPECTRA 

In this Appendix, we derive the relations between 
the three-dimensional spectrum E(k) and the one- 
dimensional spectra E^(k), EjjQc), ^ip(fc) , m ost of which 
can also be found in textbooks like Refs. [Tall7l ll8|. 



1. Total one-dimensional spectrum 



where dflk denotes the solid angle element in fc-space. 
E(k) satisfies the relation 



/ E(k)dk = 
Jo 



(A.7) 



If u is statistically isotropic in the sense that the en- 
semble average of the spectral energy of the velocity 
^|it(fc)| 2 ) is only a function of k = \k\, then E{k) be- 
comes 



E{k) = Auk' 



i\u{k)\ 2 ) 



(A. 



To evaluate Eyd in this case, we introduce cylindrical 
coordinates (k, <p, k z ) in fe-space and write the double in- 
tegral l|A.3|) in the form 



E 1D {k z ) = 2 f 
Jo 



-in 



(Nfc)l 2 ) 

2 

- (\u(k)\ 2 ) 



2ttk dn 



kdk, (A.9) 



We derive the relation between the three-dimensional 
spectrum E(k) and the total one-dimensional spectrum 
Eiu(k) = E\,{k) + 2E^{k). Consider a periodic box of 
volume V = L x L y L z with a turbulent velocity field u(x), 
which has the Fourier transform 

m= , 1 (A.i) 

with the inversion 



u(x) 



V 



j e- lk x u{k)dk 3 . (A.2) 



The one-dimensional kinetic energy spectrum is 

f f (\u(k)\ 2 ) 
E 1D (k z ) = 2 Xl V 1 dk x dk y , (k z > 0), 



(A.3) 

where (•) denotes an ensemble average, and k = 
(k x ,k y ,k z ). The factor 2 in Eq. I|A.3|) accounts for the 
fact that Eyd does not distinguish between positive and 
negative k z . Normalization of Em(k z ) is such that 



E m (k z )dk z 



V 



d.T J . (A.4) 



Equation 1|A.3|) can also be written as the a;y-average 
1 



E m (k z ) = 



L X Ly 



(\u{x,y,k z )\ 2 ) dxdy (A.5) 



and is for homogeneous turbulence equal to 
(\u(x,y,k z )\ 2 ) at any point (x,y). 

The three-dimensional velocity energy spectrum is 
given by 



E(k) ee [ 

J4tt 



V / k 2 dn k , (A.6) 



since k 2 = k 2 — k 2 z . Comparing with Eq. i|A.8|) . we see 
that 



(A.10) 



the inversion of which gives 



E(1)= . i «Ei = . ElD ^. (AJ1) 



dk 



dlnfc 



2. Longitudinal and transversal one-dimensional 
spectra 

In this section we derive the relation between longi- 
tudinal and transversal one-dimensional energy spectra 
E-^{k), and the three-dimensional energy spec- 

trum E(k). 

For homogeneous, isotropic turbulence, the energy 
spectrum tensor is given by [see, e.g.. llMll7lll8l | 



F pq (k) ee (u p (k)u*(k)) 



= [F h (k)-F T (k)\ -j^- + F T (k) S pq . (A. 12) 

If we assume incompressibility, the longitudinal compo- 
nent i<L vanishes, and thus 

F Pq (k)=(s pg -^)F T (k)=(5 pq 



and in particular 



k 2 J Airk 2 ' 
(A.13) 



**F„{k) = fl-flS, (A.14) 



k 2 ) 2k 2 

2' 



2n[F xx (k)+F yy (k)} = + (A15) 



8 



The longitudinal one-dimensional spectrum 



and thus 



£l i /.■:.) -2 // ^^dk x dk v , (A. 16) 



thus becomes 



oo oo 

E L {k z ) = J F zz (k) 2™d K =i J (l 



kl\ E(k) 



dfc, 



k 2 J k 

o k z 

(A-17) 

using the same substitution as in Eq. (|A. 9|) above. Simi- 
larly, we can write the transversal one-dimensional spec- 
trum 



E T (k z ) EE 2 JJ 



in the form 



E T (k z ) 



Taking the derivative of Eq. I|A.17|) , we find 



fc :i 



■dfc, 



(A.18) 



(A.19) 



(A.20) 



E{k) = fc 2 K( fc ) - kEi(k). 



(A.21) 



Inserting this relation into Eq. (|A.19(I allows us to express 
Er(k) through E^k) as 



E T {k) = -k 



EL(k) , £ L (fc) 



(A.22) 
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